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ABSTRACT 

A comprehensive approach is provided to the study of both S-type and P-type 
habitabihty in stellar binary systems, which in principle can also be expanded to 
systems of higher order. P-type orbits occur when the planet orbits both binary 
components, whereas in case of S-type orbits the planet orbits only one of the 
binary components with the second component considered a perturbator. The 
selected approach encapsulates a variety of different aspects, which include: (1) 
The consideration of a joint constraint including orbital stability and a habitable 
environment for a putative system planet through the stellar radiative energy 
fluxes ("radiative habitable zone"; RHZ) needs to be met. (2) The treatment of 
conservative, general and extended zones of habitability for the various systems 
as defined for the Solar System and beyond. (3) The providing of a combined for- 
malism for the assessment of both S-type and P-type habitability; in particular, 
mathematical criteria are presented for which kind of system S-type and P-type 
habitability is realized. (4) Applications of the attained theoretical approach to 
standard (theoretical) main-sequence stars. In principle, five different cases of 
habitability are identified, which are: S-type and P-type habitability provided 
by the full extent of the RHZs; habitabihty, where the RHZs are truncated by 
the additional constraint of planetary orbital stabihty (referred to as ST and PT- 
type, respectively); and cases of no habitability at all. Regarding the treatment 
of planetary orbital stability, we utilize the formulae of Holman & Wiegert (1999) 
[AJ 117, 621], as also used by previous studies, although other methods for the 
assessment of stability / instability transitions in planetary systems as, e.g., the 
Lyapunov exponent criterion have also become available. In this work we focus 
on binary systems in circular orbits. Future applications will also consider binary 
systems in elliptical orbits and provide thorough comparisons to other methods 
and results given in the literature. 
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Subject headings: astrobiology — binaries: general — celestial mechanics — 
planetary systems 



Introduction 



Starting more than a decade ago considerable observational evidence has been obtained 
indicating that planets are able to exist in stellar binary (and higher order) systems; see 



results and discussions by, e.g., Patience et al. (2002), Eggenberger et al. (2004), and Eggen- 



berger et al. (2007). These observations are in line with the empirical finding that binary (and 



higher order) systems occur in high frequency in the local Galactic neighborhood (Duquen- 
noy fc Mayor |199T Lada 2006 Raghavan et al. 2006 [Bonavita fc Desidera 2007 [Raghavan 



et al. 2010). For example, Raghavan et al. (2010) presented results of a detailed analysis of 



companions to solar-type stars, based on a sample size of 454, and concluded that the overall 
fractions of double and triple systems are about 33% and 8%, respectively, if all confirmed 
stellar and brown dwarf companions are accounted for. 

The fairly frequent occurrence of planets in binary systems is furthermore consistent 
with the presence of debris disks in a considerably large number of main-sequence star binary 
systems (e.g.. Trilling et al. 2007). In principle, as discussed by Ferryman (2011), planets 



in binary systems can be identified through two different venues: First, binaries or multiple 
star systems can be surveyed for the presence of planets by utilizing the established detection 
methods. Second, stars with detected planets can be scrutinized afterward to check if they 
possess one or more widely separated stellar companion(s); in this case, the planet (s) will 
also be categorized as belonging to a binary (or higher order) system. 

From the view point of orbital mechanics, there are two different kinds of possible orbits 
(notwithstanding positions near the Lagrangian points L4 and L5) for planets in binary sys- 



tems: S-type and P-type orbits (Dvorak 1982 ). A P-type orbit is given when the planet orbits 



both binary components, whereas in case of an S-type orbit the planet orbits only one of 



the binary components with the second component behaving as a perturbator. Eggenberger 



et al. (2004) presented a list of 15 planet-bearing binary systems with all planets in S-type 
orbits. They constitute mostly wide binaries with separation distances of up to ~6400 AU; 
however, smaller separation distances on the order of 20 AU or less have also been identified. 
In the meantime, systems with planets in P-type orbits have also been identified. Arguably, 



the most prominent case is Kepler-16, as reported by Doyle et al. (2011) and previously sug- 



gested by Slawson et al. (2011), containing a Saturnian mass circumbinary planet. Quarles 



et al. (2012) have subsequently studied this system regarding the possibility of habitable ex- 



oplanets and habitable exomoons. Recently, a transiting circumbinary multiplanet system. 
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i.e., Kepler-47, has also been identified (Orosz et al. 2012) 



There is a significant body of hterature devoted to the study of habitability in binary 
systems as well as in multiplanetary systems, which often also encompass stellar evolutionary 



considerations. Examples include the work by Jones et al. 


(200] 


), Noble et al. 


(2002 


), 


Menou 


& Tabachnik 


(2003), Asghari et al. 


(2004 


) , Sandor et al. ( i 


!007), 


Takeda et al. 


2008 




Dvorak 


et al. (2010), 


Jones & Sleep (2010) 


, and 


Kopparapu & Barnes 


(2010). An important aspect 



that has received increased recognition in the literature is that in order for habitability to 
exist and to be maintained, a joint constraint that includes both orbital stability and a 
habitable environment for a system planet through the stellar radiative energy fluxes needs 
to be met. In the framework of this paper, the zone related to this latter requirement will 
subsequently be referred to as radiative habitable zone (RHZ), which constitutes a necessary, 
though often insufficient, condition for the existence of circumstellar habitability. 

Previous work, mostly concentrated on the existence of habitability in single star multi- 
planetary systems, rendered the publication of detailed "stabihty catalogs" for the habit- 



ability zones of extrasolar planetary systems (e.g., Menou & Tabachnik 2003 Sandor et al. 
20071 iTakeda et all 120081 iDvorak et al.l 120101 iKopparapu & Barnesl 120101). For example. 



Menou & Tabachnik (2003) quantified the dynamical habitability of 85 planetary systems 
by considering the perturbing infiuence of giant planets beyond the traditional Hill sphere for 
close encounters with the theoretical terrestrial planets. They concluded that a significant 
fraction of the identified extrasolar planetary systems are unable to harbor habitable terres- 
trial planets. A statistical study on the stability of Earth-mass planets orbiting solar-mass 
stars in presence of stellar companions focusing on both the statistical properties of ejection 



times and the general prospects of planetary habitability was given by Fatuzzo et al. (2006). 



Additional work providing stability assessments for various observed extrasolar planetary 



systems based on detailed stability maps was given by Sandor et al. (2007). Takeda et 



al. (2008) explored the orbital stability or planets in double-planet systems for binaries by 
supplying an analytic framework based on secular perturbation theory; they also provided 
dynamical classification categories. Additional stability analyses to assess the habitability 



of planetary systems based on detailed numerical simulations were given by Dvorak et al. 



(2010) and Kopparapu & Barnes (2010); note that the study of Dvorak et al. (2010) also 
dealt with a limited cases of planets in double star systems in orbit either around one stellar 
component (S-type) or around both components (S-type). 

The study of planetary dynamics and habitable planet formation has meanwhile been 



described by, e.g., Quintana & Lissauer (2010) and Haghighipour et al. (2010). They show 



that Earth-mass planets are, in principle, able to form in stellar binary systems, also many 
of the physical, chemical, and astrobiological details are not fully understood, including 
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relationships between the formed planets and the orbital and physical properties of their 
host stars. The overarching conclusion of those investigations is that habitable planets in 
stellar binary (and, as anticipated, in higher-order systems) are, in general, possible, which 
is a stark motivation for providing a comprehensive study of S-type and P-type habitability 
in binary systems. The approach adopted in this study will be entirely analytic. Specifically, 
it will consider both S-type and P-type habitable orbits in the view of the joint constraint 
including orbital stability and a habitable environment for a system planet through the 
appropriate amount of the stellar radiative energy fluxes. In an earlier study, |Eggl et "ah 



(2012) focused on S-type habitability in binary systems taking into account both circular 
and elliptical orbits for the stellar binary components; this latter aspect is however beyond 
the scope of the present work as we solely focus on systems in circular orbits. Numerical 
studies for P-type habitable environments with applications to Kepler-16, Kepler-34, Kepler- 



35 and Kepler-47 have given by Kane & Hinkel (2013). 



Our paper is structured as follows: In Sect. 2, we comment on the adopted main-sequence 
star parameters and single star habitability. In Sect. 3, we introduce our theoretical approach 
suitable for stellar systems of the order of A^, although our focus will be on binary systems. 
In this regard both S-type and P-type orbits will be examined, and detailed mathematical 
criteria for the existence of S-type and P-type RHZs will be derived. In Sect. 4, we consider 
the joint constraint of planetary orbital stability and the stellar radiative environment (i.e., 
RHZs) for the establishment of circumstellar habitability. Applications regarding S-type and 
P-type systems as well as the introduction of the habitability classifications S, P, ST, and 
PT are given in Sect. 5. In Sect. 6, we convey our summary and conclusions. 



2. Stellar Parameters and Single Star Habitability 

In this study, S-type and P-type habitability is investigated mostly pertaining to stan- 
dard (i.e., theoretical) main-sequence stars. The adopted stellar parameters, which are the 
stellar effective temperatures Tes, the stellar radii i?* (which together allow to define the 



stellar luminosities L^:), and the stellar masses M^: are mainly based on the work by Gray 



(2005) (see his Table B.l) that assumes detailed photospheric spectral analyses. If a mass 
value in Gray's table appeared to be inaccurate due to numerical truncations, it was re- 
calculated if more accurate values for the surface gravity log g^: and the stellar radii were 
given. For stellar spectral types with no data available, the missing data were computed by 
employing biparabolic interpolation. 



The exception, however, are data for stars of spectral type K5 V and below. In this 
regard we relied on the results from the spectral models of R. L. Kurucz and collaborators. 
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They took into account hundreds of miUions of spectral hues for a large set of atoms and 



molecules; see Castelli & Kurucz (2004) and Kurucz (2005) for details. The effective temper- 



atures implied by these models are in close similarity to those given by Gray ( 2005 ) for most 



types of stars; however, Gray (2005) reports consistently higher effective temperatures for 



stars of spectral type late-K and M; for the latter, the difference amounts to nearly 300 K. 
Table 1 depicts the stellar parameters adopted for the present work. 

An alternative approach expected to provide very similar results for either the stellar 
luminosity or the stellar mass (with the other parameter taken as fixed) is the employment 



of a mass-luminosity relationships applicable to main-sequence stars. The work by Reid 



(1987), as well as data from subsequent studies, yield 

- = Am; 

with 7] = 0.23 and a 



2.3 for M, < 0.43 Mq and rj 



1 and a = 4.0 for M, > 0.43 Mr.-,. 



At the high- mass end, this relationship holds until about M^, = 2 Mq. It also becomes 
increasingly inaccurate for low-mass M dwarfs. Fortunately, the domains of applicability for 
Eq. (1) is consistent with most studies of binary habitability; see, e.g., Eggl et al. (2012) as 
example. 

Next we focus on single star habitability, i.e., the evaluation of various limits of habitable 



zones (HZs), which in the solar case shall be referred to as S£. Previous work by, e.g., Kasting 



et al. (1993) distinguished between conservative (GHZ) and the generalized habitable zone 
(GHZ), which can also be evaluated for general main-sequence stars, and other types of 
stars as well. For the Sun, the limits of the GHZ are given as 0.95 and 1.37 AU {i = 2 and 
4, respectively), whereas for the GHZ, they are given as 0.84 and 1.67 AU {i = 1 and 5, 
respectively); see Table 2. 



The physical significance of the various kinds of HZs obtained by Kasting et al. ( 1993 ) 



can be summarized as follows: The GHZ is defined as bordered by the runaway greenhouse 
effect (inner limit) and the maximum greenhouse effect (outer limit). Goncerning the latter 
it is assumed that a cloud-free GO2 atmosphere is still able to provide a surface temperature 
of 273 K. The inner limit of the GHZ is defined by the onset of water loss. In this case, a 
wet stratosphere is assumed to exist where water is lost by photodissociation and subsequent 
hydrogen escape to space. Furthermore, the outer limit of the GHZ is defined by the first 
GO2 condensation attained by the onset of formation of GO2 clouds at a temperature of 



273 K; see, e.g.. Underwood et al. (2003) and Selsis et al. (2007) for further details. Table 3 



conveys the results for the HZs of the different type of main-sequence stars considered in the 
present study. 



For the outer edge of circumstellar habitability, even less stringent limits have been 



introduced in the meantime (e.g., Forget & Pierrehumbert 1997 Mischna et al. 2000). They 
are based on the assumption of relatively thick planetary CO2 atmospheres as well as strong 
backwarming that may further be enhanced by CO2 crystals and clouds. These limits, which 
in case of the Sun correspond to 2.4 AU {se = sq), conform to the extended habitable zone 
(EHZ), have also been taken into account in our study, although the significance of the EHZ 
has meanwhile been criticized as a result of detailed planetary radiative transfer models 
(Halevy et al. 2009). Moreover, in the framework of the present study, we also consider 

S3; 



planetary Earth-equivalent positions defined as -Re,eqv — y-^*/-^© and labelled as se 
see Table 3. It is meant as an intriguing reference distance of habitability both regarding 
single stars and stellar binary systems. 



3. Theoretical Approach 

3.1. Basic Equations 

Next we introduce the governing equations for investigating the RHZs of binary systems 
pertaining to both S-type and P-type orbits. This approach targets the requirement of 
providing a habitable environment for a system planet based on the radiative energy fluxes 
of the stellar components. The requirement of planetary orbital stability will be disregarded 
for now; it will be revisited in Sect. 4. The importance of orbital stability for allowing 
circumstellar habitability in stellar binaries will, however, be considered in an appropriate 
and consistent manner in the main body of the study. 

For a star of luminosity Li, given in units of solar luminosity Lq, the distance di of the 
habitability limit si as identified for the Sun, which may constitute either an inner or outer 
limit of habitability (except £ = 3), is given as 



di = se^l ^\ (2) 

*Jrel,i-'^0 

In case of a multiple star system of order N with distances di, the limit of habitability related 
to Se is given a^ 

^ Li Lq 

i=l "-^reljiU-i i'^ 



^ This equation is analogous to an equivalent equation of electrostatics relating a general distribution of 
charges to the resulting electrostatic potential in free space (Jackson 1999 see p. 40, Eq. (1.48)); a modified 
version of Eq. (3) has previously been considered by, e.g., Eggl et al. (2012 1. 
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In Eq. (2) and (3), Sj.c\,i = 5'rci,j(Tcfj) describes the stellar flux in units of the solar 
constant that is a function of the stellar effective temperature T^s (e.g., Kasting et al. 1993 



Underwood et al. 2003). Specifically, using the formalisnJ by Selsis et al. (2007), we find 



that 



rel,i 



(4) 



with S£, a^, and bz in AU, and = T^s — 5700 in K. Selsis et al. (2007) also found that 



for Se < 1, corresponding to inner limits of habitability, the fitting parameters are given as 
Qz = 2.7619 X 10~^ and bz = 3.8095 x 10~^, whereas for sg > 1, corresponding to outer limits 
of habitability, they are given as = 1.3786 x 10""^ and bz = 1.4286 x 10~^; note that se = 1 
corresponds to the customary notion of Earth-equivalent positions. Appropriate values for 
Si are given in Table 2. 

In the following, we will focus on the case of binary systems, i.e., N = 2. In this case 
Eq. (3) reads 

■ ^ (5) 



5'rel,l'^l 



+ 



5'rel,2C^2 



with 



dl 



+ z'^ + 2az cos (p 
+ z'^ — 2az cos ip 



(6a) 
(6b) 



Here a denotes the semidistance of binary separation, z the distance of a position at the 
habitability limit contour (which later on will be referred to as "radiative habitable limit", 
see below), and if the associated angle; see Fig. 1 for information on the coordinate set-up 
for both S-type and P-type orbits. We will also assume Li > L2 without loss of generality. 



With L'- defined as 



L' 



0^rel,i 



(7) 



^Here S'rci.i represents the normalized stellar flux in units of the solar constant, 1368 W ni~^, given by the 
stellar spectral energy distribution. Therefore, ordinarily, no si dependence for S'rei,! should exist. However, 
the formulae by Selsis et al. ( 2007 ) utilize previous results by Kasting et al. ( 1993 ) who provided numerical 
values for limits of habitability for different types of stars considering various limit definitions (i.e., se values 
identified for the Sun). But Kasting et al. (19931 used for the solar effective temperature an unusually low 
value of 5700 K instead of 5777 K as currently accepted (e.g., Stix|2004 ). Hence transforming the polynomial 
fit based on the work by Kasting et al. aimed at considering the correct solar effective temperature renders 
a weak dependence on for the Sre\,i values. In contrast the method by Underwood et al. ( 2003 ) provides a 
polynomial fit for Srei,i without considering the solar temperature revision. An alternative method has been 



used by Cuntz & Yeager ( 2009 ) and subsequent work. In this approach the polynomial fit by Underwood 



et al. ( 2003 ) is corrected via a triangular function based on data for stars of spectral type FO V, GO V, and 
KO V. As a result the corresponding S'roi,i values do also not show any si dependence. 
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henceforth referred to as recast stellar luminosity, z{ip) is given as 

+ A2Z^ + Aiz + Ao = 



with 



A2 = 2a2(l - 2 cos^ if) - sj{L[ + L'^) (9a) 
Ai = 2as1 cos Lp{L\ — L'2) (9b) 
Aq = a^- a^sl{L\ + L'^) . (9c) 



Equation (8) constitutes a fourth-order algebraic equation that is known to possess four 



possible solutions (Bronshtein & Semendyayev 1997), although some (or all) of them may 



constitute unphysical solutions, i.e., z{lp) having a complex or imaginary value. The adopted 
coordinate system constitutes, in essence, a polar coordinate system except that negative 
values for z are permitted; in this case the position of z is found on the opposite side of angle 

In principle, it is possible to consider for Eq. (8) to only have solutions for z{ip) given 
as z > 0; in this case the entire interval for z{ip), which is < < 27r, needs to be 
examined. The following types of solutions are identified: For S-type orbits, two solutions 
exist in the intervals centered at = and at = vr (or one coinciding solution at each 
tangential point); see Fig. 1. However, there will be no solution in the typically relatively 
large intervals containing (p = ti/2 and (p = 3tt/2. Clearly, the size of any of those intervals 
critically depends on the system parameters a, L[, L'2, and si, as expected. For P-type orbits, 
on the other hand, there will be one solution for each value of (f in the range of < < 27r. 

However, in general, negative values for the solutions of z{ip) also exist. If taken into 
account, it will be sufficient to restrict the evaluation of Eq. (8) to the range < ^9 < vr. In 
this case, for S-type orbits, there will be four solutions in the intervals with endpoints ip = 
and 99 = vr, as well as two solutions (if L[ 7^ Lg) in a more extended interval containing 
these points. Also, a pair of solutions will become one coinciding solution at each tangential 
point. However, again, there will be no solution in the interval containing ip = 7t/2. In case 
of P-type orbits, there will be two solutions for any value of p in the range of < p < 71. 
We will revisit this assessment in conjunction with the algebraic method for attaining the 
solution zf, additionally, detailed mathematical criteria will be given for the existence of 
RHZ for S-type and P-type orbits. 



Next we will focus on equal-star binary systems. Detailed solutions for general binary 
systems (i.e., systems of stellar components with by default unequal masses, luminosities. 
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and effective temperatures) pertaining to both S-type and P-type orbits will be given in 
Sect. 3.3. Both subsections will be aimed at deriving RHZs; see, e.g., [Williams fc Pollard" 



(2002) for general discussions on the role of RHZs for the attainment of habitability in 
star-planet systems. However, strictly speaking, they will deal with identifying radiative 
habitable limits (RHLs) connected to a distinct value of noting that manifesting a RHZ 
requires that the RHL for s^^out to be located completely outside of the RHL for s^^in with 
se^iri and s^^out appropriately paired. A summary about the existence and structure of the 
RHZs, encompassing the radiative CHZs, GHZs, and EHZs, will be given in Sect. 3.4; this 
subsection will also convey cases where no RHZs exist due to the behavior of the RHLs 
regarding se^in and s^,out- 

However, in none of these subsections the requirement of orbital stability for possible 
system planets is considered, which in general entails further restrictions on the manifestation 
of binary system habitability. This latter aspect will be discussed in Sect. 4 and considered 
in detail in the subsequent set of applications (see Sect. 5). 



3.2. Equal-Star Binary Systems 

Now we focus on the special case of equal-star binary systems, i.e., stars of identical 
recast luminosities, i.e., L'^ = L'2 = L'. For theoretical main-sequence stars this assumption 
also implies S'rei,i = 5'rei,2 and Mi = M2; this latter assumption about the stellar masses is 
relevant for the orbital stability constraint of system planets. With Ai = 0, Eq. (8) now 
constitutes a biquadratic equation that can be solved in a straightforward manner. The 
other coefficients are given as 

A2 = 2a^{l - 2 cos^ ip) - 2s]L' (10a) 
Aq = a^- 2a^sjL' . (10b) 

Thus, the solution of Eq. (8) is given as 



z = ±\/-a^{l-2cos'^ip) + sjL'±^/D2 (11) 

with 

D2 = 4L'' + 4a^cosV-4a2cosV(a2-s2L'). (12) 

With known systems parameters, which are a, se, and L', the function z{(p), describing the 
habitability limits for the binary system associated with inner limit and outer limit values 
Se derived for the Sun (see Sect. 2 for details) can be obtained in a straightforward manner. 
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Owing to the system symmetry, the existence of S-type and P-type RHLs can be found 
by attaining the solutions of Eq. (11) for = and </? = 7r/2. First we examine the solutions 
of Eq. (11) for (fi — 0, i.e., cost/? = 1, which are given as 

z = ±^ a'^ + sjU ± se^ s'jL'^ + Aa'^L' . (13) 

The total number of solutions for Zi is four as expected; it collapses to two solutions if the 
value of the interior square root approaches zero, which obviously reflects the mathematical 
hmit for the existence of S-type RHLs. Thus, S-type RHLs exist if 

a > seV2IJ . (14) 

Next we examine the solutions of Eq. (11) for = 7i/2, i.e., cosip = 0. This allows us to 
gauge the existence of P-type RHLs; the latter implies two solutions of Eq. (11) regardless 
of the value for (p. Here the solution is given as 

z = ±^-a? + 2s}U . (15) 
Thus, the condition for P-type RHLs is given as 

a < Si\/2L' . (16) 

Therefore, Eqs. (14) and (16) allow to identify conditions for S-type and P-type RHLs, 
respectively, for equal-star binary systems, which depend on the systems parameters a, s^, 
and V . Comparing Eqs. (14) and (16) also implies that the joint existence S-type and P-type 
habitability in equal-star binary systems is impossible, irrespectively of the planetary orbital 
stability requirement (see Sect. 4), noting that the latter imposes an additional constraint 
on habitability even when the RHZ- related conditions are met. Figure 2 depicts the borders 
of the S-type and P-type radiative habitable limits, i.e., RHLs, for different values of in 
regard to a and L' . 
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3.3. General Binary Systems 

3.3.1. Method of Solution 
We now focus on obtaining solutions for our key equation, Eq. (8), pertaining to S-type 



and P-type RHLs for general binary systems, i.e., L'^ ^ L'2. Following, e.g., Abramowitz & 



Stegun (1972) and Beyer (1987), the set of solutions for a fourth-order polynomial reads 



zz 

Zi 



2 2 

2 2 

+-C - -8 
2 2 

^2" ^2' 



wi 



-2a2(l - 2 cos2 y?) + s}{L'-^ + L'^) + yi 



C = 

V = -^sKL; + ^'2) +4as2(L; -L'2)C-icos(/?-4a2(l -2cosV) -Z/i 
£ = ^Jsj{L[ + L'2) - 4:asj{L[ - L'2)C-^ cosip - 40^(1 - 2cos2 ip) - yi 

with yi as a solution of the resolvent cubic equation 

y' - A^y^ - Aoy + {4A2A0 - Aj) = 



(17a) 
(17b) 
(17c) 
(17d) 

(18a) 
(18b) 
(18c) 

(19) 



with Aq, Ai, and A2 given by Eqs. (9a) to (9c). The solutions Zi given through Eqs. (17a) 
to (17d), if existing, are ordered as Zi < Z2 < Z3 < for ip = 0, as well as for any other value 
of ip as subsequently identified. Although Eq. (19) has three possible solutions, there is only 
one appropriate choice for yi because it is necessary to avoid that all Zi obtained through 
Eq. (8) are of imaginary or conjugate complex value in cases where S-type or P-type RHLs 
exist. 

The adopted method also requires a solution for y; see Eq. (19). The only acceptable 
solution for y, customarily denoted as yi, is given as 



yi = -^^2 + {S + T) 



(20) 



^In Eqs. (18b) and (18c) we also correct a typo in the sign preceeding the term of C ^. 
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with the substitutions 



R+VD3 (21a) 



and with Q and R given as 



T = yn-^/Ds (21b) 

L>3 = g^ + i?' (21c) 



Q = ^«i-^«2 (22a) 
R = -^oo + ^oias - (22b) 



while noting that 



do = - Al (23a) 

61 = -4Ao (23b) 

62 = -A2 . (23c) 

These sets of equation can be solved and appropriate values for z{(p) can be obtained. The 
results will depend on the system parameters a, L[, L2, and se, as expected. 

Next we describe the solutions for S-type and P-type RHLs in more detail. It is impor- 
tant to recognize that apriori choices about the existence of S-type and P-type RHLs are 
neither necessary nor possible as the existence of any of those RHLs is determined by the 
fulfillment of well-defined mathematical conditions; they will also be given in the following. 



3.3.2. S-type Orbits 

An analysis of the possible solutions for Eqs. (17a) to (17d) shows that for S-type RHLs 
valid solutions are obtained based on 

S = ^/R^ + K^ (cos^ + isinO (24a) 
T = \/R^ + K'^ (cos^-isinO (24b) 



with 



K = ^/\D^\ (25a) 
e = ^arctan(^) (25b) 
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with R given by Eq. (22b). Therefore, the solution of the resolvent cubic equation, Eq. 
(19), is given as 

yi = -\a2 + 2\^Wn^cosC (26) 

For the values of Zi it is found that for the interval centered aX = 0, Zi and Z2 exhibit 
negative values, whereas Zo, and z^ exhibit positive values with Zi < Z2 < z^ < z\. Thus, z\ 
and Z2 describe the RHL regarding star SI, whereas Zj, and z^ describe the RHL regarding 
star S2 (see Fig. 1). Conversely, for the interval centered at 9? = tt, Z\ and Z2 again exhibit 
negative values and zj, and z^ exhibit positive values. In this case, z\ and Z2 describe the 
RHL for star S2, whereas Zz and 2^4 describe the RHL for star SI. Surely, no solutions are 
obtained in the vicinities oi — 7r/2 and </? = 37r/2, as expected. Thus, for each angle </? in 
the range of < </9 < 27r the appropriate number of solutions is attained to describe S-type 
RHLs. However, due to symmetry, solutions are only needed for < </? < tt. 

The existence of S-type RHLs requires that Z2 < Z3 because otherwise the two distinct 
S-type RHL contours about the two binary components would not be separated, which 
corresponds to the condition 

C >\{V + £). (27) 

This equation can most likely be appropriately rewritten to provide an expression only 
containing the system parameters Aq, Ai, and A2 defined through Eqs. (9a) to (9c). However, 
the resulting expression is expected to be highly intricate. 



3.3.3. P-type Orbits 

An analysis of the possible solutions for Eqs. (17a) to (17d) also shows that for P-type 
RHLs valid solutions require 

L>3 > ; (28) 

see Eq. (21c). The detailed evaluation of this condition requires the evaluation of various 
sets of equations denoted as Eqs. (9), (21), and (22); see Sect. 3.1 and 3.3.1. 

In terms of the solutions for P-type RHLs it is found that for < (/? < 71/ 2 and 
37r/2 < ip < 271, zi exhibit negative values and Z2 exhibit positive values, whereas z^ and Z4 
are undefined; they are also not needed for outlining P-type RHLs. Moreover, for the range 
of 7r/2 < (p < 37r/2, z^ exhibit negative values and Z4 exhibit positive values, noting that 
zi and Z2 remain undefined. For 7r/2 and 37r/2, removable singularities are identified, which 
can easily be fixed through interpolation taking values of z for neighboring angles of (f. In 
summary, for each angle (p in the range of < </? < 27r two values of Zi (i.e., one positive 
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and one negative value) are identified allowing to determine P-type RHLs. However, due to 
symmetry, solutions are only needed for < (/? < vr. 

Furthermore, through analytical transformations it can be shown that the condition 
depicted as Eq. (29) can be rewritten as 

IGAoAj + U4:AoAlA2 - USAlAj + 256AI - AAjA^ - 11 A\ < . (29) 

with Ao, Ai, and A2 defined through Eqs. (9a) to (9c). A key aspect concerns the following: 
For P-type RHLs to exist for the system parameters a, L[, L2, and se, it is necessary that 
the relations (28) and (29), which are equivalent, must be fulfilled for any angle of (p. 

A subsequent evaluation of Eq. (29) shows that in case of equal-star binary systems, 
corresponding to Ai = 0, Eq. (29) can be simplified as 

(A2-4Ao)'>0; (30) 

this relationship is fulfilled in a trivial manner. However, if the positive root of Eq. (30) is 
taken, the then recovered relation between Aq and A2 leads to the relationship of Eq. (16) 
valid for the special case of equal-star binary systems previously discussed (see Sect. 3.2). 



3.4. Calculation of RHZs for Binary Systems 



The identification of the RHZs in binary systems requires the calculation of limits of 
habitability, i.e., RHLs, as discussed in Sect. 3.1. They need to be established for values of si 
with i = 1, 2, 4, 5, and 6 (see Sect. 2 and Table 2), which are informed by model-dependent 



physical limits of habitability for the solar environment (e.g., Kasting et al. 1993). As part 



of the process, the parameters of se need to be appropriately paired in terms of the inner and 
outer limits of habitability. For the CHZ the parameters (s^^in, s^^out) need to be paired as 
(s2, S4), whereas for the GHZ, they need to paired as (si, S5). For the EHZ, the parameters 
of Se need to paired as (si, Sq), considering that both the CHZ and the GHZ shall be viewed 
as subdomains of the EHZ. 

For S-type and P-type orbits, the radiative zones of habitability IUiZ{z), which consti- 
tutes a circular region (annulus) around each star SI and S2 (S-type) or both stars (P-type) 
can be determined as 



RHZ(z) = Mm(n{z,a) 



Max(n{z,a] 



and 



Rm{z) = Mmfniz,^)) -Max(7^(2,^) 



(31) 
(32) 
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respectively; see Fig. 1 for coordinate information. Here TZ{z, a) and TZ{z, (p) describe the 
areas bordered by the RHLs defined by s^ in and s^^out- The calculation of the extrema is 
applied to the angles a and ip for the intervals < a < vr and < ip < 7i/2, respectively; 
note that we assumed Li > L2 without loss of generality. In the S-type case the calculation 
of the extrema pertaining to the RHZ values is based on the angular coordinate a instead 
of ip; however, the angular coordinate (p is still needed for the calculation of Zi as part of the 
overall approach toward identifying S-type and P-type habitability. 

Figure 3 depicts examples of RHLs and RHZs for different types of systems. In the S- 
type case the RHLs arc bended toward the center of the system, whereas in the P-type case 
they arc of notable elliptical shape. The RHZs always constitute circular annuli obtained 
through inspecting the appropriate minima and maxima of the RHLs. The examples as 
depicted include S-type and P-type systems with separation distances (2a) of 0.5 AU and 
5.0 AU, respectively. Cases of both equal-star and non-equal star binaries are selected. The 
focus of this figure is the identification of the appropriate circular region (annulus) for each 
case. The figure also indicates the portions within the TZ{z, a) and TZ{z, (p) domains that are 
not part of the RHZ (2;) annuli. 

Next we determine the values for the extrema pertaining to RHZ(z) following Eqs. (31) 
and (32) based on the solutions of Eq. (8) given as Eqs. (17a) to (17d). Noting that for 

existing values of Zi they are ordered as Zi < Z2 < z^ < Z4^, and negative values for Zi are 
permitted it is sufficient for S-type orbits, both for the stars SI and S2, to only consider 
solutions for (/? = 0. For P-type orbits, a more detailed assessment is required (see below). 
For S-type orbits, regarding star SI, the extrema are obtained as 



se,out 



|a + 22(0)1 
la + zi 



RHZin = Max{n{z,a)) 
RHZout = Mm{n{z,a)) 
and for star S2, they are obtained as 

RHZin = Max{n{z,a)) 

RHZout = Mm{n{z,a)) 

The size of each annulus Aj^jjz fo^ pairs (s^^in, -s^^out) is given as Arhz — RHZout — RHZ 
For P-type orbits, the extrema are given as follows 



\a - Z3 
|a- 24(0)1 



St out 



(33a) 
(33b) 

(34a) 
(34b) 



RHZin = Max{n{z,cp)) 
RHZout - Mm{n{z,ip)) 



= l^i(0)| 



|^l(<^OUt)| 



(35a) 
(35b) 
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We note that for the angle (pont for R-HZout no straightforward expressioiij^ exists; it is located 
in the interval < ipout < tt/2. But it can be found numerically as it is given by the angle 
where the minimum of |-z(v5)| occurs. 

There is also another complication in the identification of RHZ P-type orbits. Generally, 
it is required for the RHL of s^ out to be located completely outside of the RHL of s^^in, 



Mm(n{z,ip)) > Max(n{z,(p) 

^ ^ out ^ 



I.e., 



(36) 



However, this condition is violated in some models, especially for relatively large values of a 
as well as relatively small ratios of L'2/L'^. In this case the RHZ for (s^^in, s^^out) is nullified, 
a behavior that may occur for the pairings (52,^4), (^1,55), and (si,S6), corresponding to 
the GHZ, GHZ, and EHZ, respectively. In this regard the existence of the GHZ is in most 
jeopardy as (S4 — S2) constitutes the smallest bracket among the various kinds of HZs (see 
Table 2). Also note that if the GHZ is nullified, the GHZ will be nullified as well considering 
that the GHZ (if existing) is entirely located within the GHZ. Likewise, if the EHZ is nullified 
the existence of both the GHZ and the GHZ will be nullified. Detailed examples will be 
given in the application segment of this paper; see Sect. 5.2 and 5.3 for details. However, 
this phenomenon does not exist for RHZ S-type orbits. The reason is that in comparable 
settings the RHL for s^^out is mathematically undefined, which also amounts to a limitation 
for the extent of the associated RHZ. 

For equal-star binary systems, with the property of L'^ = Lg, expressions for RHZin 
and RHZout for S-type and P-type orbits can be obtained based on Eqs. (13) and (15). For 
S-type orbits we find 



RHZir 



RHZ, 



+ 4out^' - Stout a/ sLut^'' + 4a2L 



'^.out-" 



(37a) 
(37b) 



It is furthermore intriguing to explore the limit a ^ S£_outV^, noting that s^^outV^ > 
se iny/L'- In this case it is found that 



RHZin 

RHZout 



2a 
1 

2a 



'out 



(38a) 
(38b) 



^Due to the nature of the underlying equations for Zi, an analytic expression for (p is deemed possible. 
However, it is expected to be highly complicated and thus left as a student exercise. 
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Moreover, in the limit of a — )■ oo, the expressions for single star HZs regarding RHZjn and 
RHZout are recovered, as expected, which are given as s^^^^^L' and s^^outV^, respectively. 



They are in agreement with the expressions previously obtained by Kasting et al. (1993) 



Underwood et al. (2003), Selsis et al. (2007), and others. 



Results for P-type orbits can be obtained considering 



RHZin 

RHZout 
In this case we find 



Max(^7^(2,0),7^(^ 
Minf7^(z,0),7^(z, 



TT, 

2- 

-) 



7^(2,o) 



RHZi 



0? + s]-L' + Sr, 



s2 „L'2 + 4a2L' 



RHZ 



out 



(39a) 
(39b) 

(40a) 
(40b) 



based on the system parameters a, L', s^^in, and s^^out- Additionally, the requirement to 
avoid that the RHL for s^^out to be partially or completely located inside of the RHL for s^^in, 
see Eq. (36), entails 



a2 + s,,inWa2L'+-4,,L'2 < 



+ r,ufL 



(41) 



which allows to set constraints on the separation distance 2a of the binary system noting that 
the values of L', s^^m, and s^^out are all subject to distinct restrictions, particularly in case of 
main-sequence stars (see Tables 1 and 2). Depictions of the condition (41) for the pairings 
(-52,54), ("Si,S5), and {si,Sq), corresponding to the CHZ, GHZ, and EHZ, respectively, are 
given in Fig. 4. A similar expression is expected to hold for nonequal-star binaries, though 
it will be highly complicated. 



4. Constraints on Habitability due to Planetary Orbital Stability 

A primary constraint on planetary habitability is that planets are required to exist in 
the HZ for a sufficient amount of time allowing basic forms of life to emerge and develop. In 
order to adhere to this criterion, planetary orbital stability is required. There is a significant 
body of hterature devoted to this topic, including studies of binary and multi-planetary 



systems, which often also consider aspects of stellar evolution (e.g., Jones et al. 


2001 Noble 


et al. 2002; Menou & Tabachnik 2003 Sandor et al. 2007; 


Takeda et al. 2008 I 


Dvorak et al. 


2010 Haghighipour et al. 2010 Kopparapu & Barnes 201C 


)). 
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Early studies of planetary orbital stability pertaining to planets in both S-type and 
P-type orbits demonstrated that planets can exist in systems of binary stars for 3000 binary 



periods Dvorak (1984, 1986). Although these investigations considered relatively short in- 



tegration times, Dvorak determined upper and lower bounds of planetary orbital stability 
considering the orbital elements, semimajor axis and eccentricity, of the proposed binary 
stars. Since this pioneering work, many additional studies have been performed. The fore- 
most investigation extended the original study by a factor of 10 in integration times and an 
extended range of orbital elements ( Holman &: Wiegert||1999 ). In addition, the nature of the 
bounding formula was derived and discussed using a more statistical framework 



Holman 



&: Wiegert] developed fitting formulae for both S-type and P-type planets in binary systems 
given as 

0.464 - 0.38/1 + J's(/i,eb) (42) 



QjQJ- 

a 



and 



respectively. 



Oct 

a 



1.60 + 4.12/i + J'p(/i,eb) 



(43) 



These equations give the critical semimajor axis in units of the semimajor axis a in 
case of S-type and P-type orbits. For an S-type orbit, the ratio flcr/a, see Eq. (42), conveys 
the upper limit of planetary orbital stability, whereas for a P-type orbit, the ratio acr/a, 
see Eq. (43), conveys the lower limit of planetary orbital stability. Moreover, /i denotes the 
stellar mass ratio given as /i = M2/(Mi+M2), where Mi and M2 constitute the two masses of 
the binary components with M2 < Mi. Equations (42) and (43) also contain the parameter 
functions J-s(/i, Cb) and J^p(/i, Cb), which depend on the aforementioned mass ratio /i and the 
eccentricity of the stellar binary, Cb. Considering that this paper is solely aimed at stellar 
binaries in circular orbits (i.e., Cb = 0), it is found that J-'s = J-'p = 0. 

Planetary orbital stability has been investigated by many authors using chaos indicators, 
such as the maximal Lyapunov exponents (MLE), fast Lyapunov indicator (FLI), and the 
mean exponential growth factor of nearby orbits (MEGNO), to name those commonly used; 
see, e.g., Satyal et al. (2013) for details, recent applications, and references. These methods 



have also been used to characterize the transition from stable to unstable orbits within the 



framework of the circular and elliptical 3-body problems; see, e.g., Cuntz et al. (2007), Eberle 



et al. (2008), and Szenkovits & Mako (2008) for details. 



In our previous work (e.g., Musielak et al.|[2005 ), we studied the stability of both S-type 
and P-type orbits in stellar binary systems, and deduced orbital stability limits for planets. 
These limits were found to depend on the mass ratio between the stellar components and 
the distance ratio between planetary and binary semimajor axes. This topic was revisited 



by Eberle et al. (2008), who used the concept of Jacobi's integral and Jacobi's constant to 
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deduce stringent criteria for the stability of planetary orbits in binary systems for the special 
case of the coplanar circular restricted three-body problem. Recently we studied planetary 
orbital stability through the perspective of a chaos indicator, the MLE ( Quarles et al.|[20TT ). 
From the use of a chaos indicator a cutoff value for the maximum Lyapunov exponent was 
determined as an additional stability criterion for S-type planets in the circular restricted 
3-body problem. 



5. Applications and Classifications 

5.1. Case Studies of S-Type Habitability in Binaries 

In the following, we provide selected case studies of S-type habitability for binary sys- 
tems, including systems of equal and non-equal masses (see Table 4). Our intent is to 
demonstrate the functionality of the method-as-proposecQ there is no intent to offer de- 
tailed parameter studies of systems or give comparisons to other work; these efforts will be 
pursued in later papers of this series. 

Figure 4 allows comparative insight into S-type habitability for binary systems. The 
selected cases are systems with masses of Mi = M2 = 1.0 Mq and Mi = 1.5 Mq, M2 = 
1.0 Mq with binary separation distances of 10 AU and 20 AU. For single stars of 1.0 Mq, 
the GHZ extends from 1.049 to 1.498 AU, and the GHZ extends from 0.927 to 1.831 AU; 
these values are modestly higher than those for G2 V stars given in Table 3 owing to a 
minuscule difference in mass (i.e., 0.99 versus 1.0 Mq). For an equal-mass binary system of 
1.0 Mq with a separation distance (2a) of 10 AU, the radiative GHZ and GHZ extend from 
1.056 to 1.511 AU, and from 0.932 to 1.853 AU, respectively, for each component; the extent 
of the radiative EHZ is altered from 2.64 to 2.70 AU. However, there is now an (upper) 
orbital stability limit of 1.37 AU imposed on each star; consequently, significant portions of 
the radiative GHZ and GHZ are unavailable for providing habitability. If the second star is 
placed at a distance of 20 AU, the changes in extent for the radiative GHZ and GHZ relative 
to the single star case are extremely minor. Specifically, for binary separations of 10 AU 
and 20 AU, the sizes of the radiative GHZ increased by 1.9% and 0.6% relative to the case 
of a single star case. Moreover, for the system with a separation distance of 20 AU, the 
imposed orbital stability limit is at 2.74 AU; consequently, the full extents of the radiative 
GHZ, GHZ, and EHZ are now available for planetary habitability. 



'^The method has also successfuUy been used for determining the radiative habitable zone of Kepler-16, 
a binary system with Mi = 0.69 Mq and M2 = 0.20 Mq, and a Saturnian planet in a P-type orbit; see 
Quarles et al. (2012) for a detailed study of the system's habitability. 
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Figure 4 also shows results for the pairs Mi = 1.5 Mq and M2 = 1.0 Mq. In case of a 
single 1.5 Mq mass star, the CHZ and GHZ extend between 1.88 and 2.49 AU, and between 
1.65 and 3.11 AU, respectively, whereas the EHZ extends up to 4.61 AU. In this type of 
system, a secondary star of 1.0 Mq placed at a separation distance of 10 AU again modifies 
the extents of the radiative CHZ and the radiative GHZ, which are now given as from 1.89 
to 2.49 AU and from 1.66 to 3.14 AU, respectively; however, there is an imposed orbital 
stability limit at 1.56 AU. Therefore the entire extents of the radiative CHZ, GHZ, and EHZ 
of the primary are unavailable to the facilitation of planetary habitability. If the secondary 
star is placed at a separation distance of 20 AU, the radiative CHZ, GHZ, and EHZ of the 
primary star are again similar to those of a 1.5 Mq mass star. However, the orbital stability 
limit now occurs at a distance of 3.12 AU from the primary; therefore, the entire radiative 
EHZ (if measured from the outer edge of the GHZ) is available for habitability. 

In summary, for potentially habitable S-type binaries, due to the nature of the relatively 
large separation of the stellar components, the effect of the stellar secondary on the extents 
of the RHZs is usually very minor, i.e., about a few percent or less, with the biggest impact 
occurring in F-type systems. The secondary's main influence on circumstellar habitability 
consists in imposing restrictions on the facilitation of habitability due to its gravitational 
influence, which is most relevant in close binary systems; see Sect. 5.3 for additional results. 
The results of the present work refer to the special case of binaries in circular orbits. More 
detailed results, which also encompass binary systems in elliptical orbits, thus entailing much 



more intricate outcomes, have been given by, e.g., Eggl et al. (2012); these types of studies 
will be revisited in Paper II. 



5.2. Case Studies of P-Type Habitability in Binaries 

Various sets of models have also been pursued for P-type habitability (see Fig. 6). In 
this case, we focus on systems with masses of Mi = M2 = 1.0 Mq and Mi = 1.5 Mq and 
M2 = 0.5 Mq; additionally, we also considered Mi = 1.25 Mq and M2 = 0.75 Mq (see 
Tables 4, 5, 6, and 7 for details). The separation distances (2a) were chosen as 0.5, 1.0, and 
2.0 AU. Our approach is again a two-step process. We first explore the existence and extent 
of the radiative CHZ, GHZ, and EHZ; subsequently, we consider the additional constraint of 
planetary orbital stability, which in the case of P-type orbits constitutes a lower limit (see 
Sect. 4). The results for the occurrence of habitability can be summarized as follows. 

For systems with masses of Mi = M2 = 1.0 Mq, the following behavior occurs: For 
separation distances (2a) of 0.5 AU, the inner limit (i.e., RHL; see Sect. 3.4) of the radiative 
CHZ {si = 0.95 AU) is found to vary between 1.46 and 1.54 AU, with 1.54 AU being the 
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acceptable inner limit value; see Eq. (32). Furthermore, the outer limit of the radiative GHZ 
[sg = 1.37 AU) varies between 2.10 and 2.16 AU, where 2.10 AU constitutes the acceptable 
limit value. In consideration of the orbital stability limit at 0.915 AU, the entire extent of the 
radiative GHZ offers circumbinary habitability. The acceptable inner limit of the radiative 
GHZ is given at 1.38 AU, whereas the acceptable outer limit occurs at 2.58 AU; thus, also 
the entire radiative GHZ is identified as habitable. 

For separation distances of 1.0 AU, the orbital stability limit is given at 1.83 AU, which 
falls inside the domain of the radiative GHZ ranging from 1.69 to 2.06 AU; therefore, only 
half of the radiative GHZ is available for circumbinary habitability, whereas the other half 
is not. Noting that the radiative GHZ is fully embedded into the radiative GHZ, entails 
that also only a fraction of the radiative GHZ is available for circumbinary habitability. 
However, the full extent of the radiative EHZ, which has an inner limit at 3.70 AU and an 
outer limit at 3.83 AU is able to offer habitability, if existing. We also studied models with 
binary separations of 2.0 AU. In this case, the orbital stability limit is found at 3.66 AU. 
Therefore, both the radiative GHZ and GHZ are unavailable for providing habitability; the 
latter has an outer limit that varies between 2.39 and 3.05 AU with 2.39 AU constituting 
the appropriate value. The radiative EHZ is found to vary between 3.60 and 4.09 AU. With 
3.60 AU to be ruled acceptable as the conservatively selected (i.e., inner) limit of the EHZ, 
none of the radiative EHZ should be considered acceptable for circumbinary habitability. 

Further case studies were pursued for systems with unequal distributions of masses, 
and by implication, luminosities as, e.g., systems with Mi = 1.5 Mq and M2 = 0.5 Mq. 
Following the well-known mass-luminosity relationship for main-sequence stars, it is found 
that a 1.5 Mq star possess a luminosity L^, about 3.5 times higher than a 1.0 Mq star (see 
Table 4); a similar factor of difference applies to the recast stellar luminosity L' (see Sect. 3.1). 
Thus, the combined luminosity of the (1.5 Mq,0.5 Mq) system is considerably higher than 
the combined luminosity of the (1.0 M0,1.O Mq), as expected, and the distribution of 
luminosities is highly uneven. On the other hand, following the formula of [Holman &: Wiegert 



(1999), as well as previous work by Dvorak (1986), an unequal distribution of stellar mass 



i.e., a smaller value of /i (see Sect. 4), entails a smaller orbital stability limit. Because it 
constitutes a lower limit, i.e., positioned relatively close to the two stars, it offers a relatively 
large "window of opportunity" for planets in any of the RHZs to be orbitally stable. 

Results are given in Fig. 6. For a binary separation of 0.5 AU it is found that both the 
radiative GHZ and GHZ exist, and they are both fully permitted based on the planetary 
orbital stability constraint, although the width of the GHZ is relatively small. The GHZ 
extents from 2.14 to 2.28 AU, whereas the GHZ extents from 1.91 to 2.90 AU; the orbital 
stability limit is given at 0.66 AU. There is also a considerably large domain of the EHZ, 
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which has an inner limit at 4.40 AU and an outer limit at 4.89 AU. 

However, for binary separation of 1.0 AU, the radiative CHZ has been nullified, even 
though its existence would be consistent with the planetary orbital stability constraint con- 
veying a lower limit at 1.32 AU. The reason is that the RHL of S4, which is at 2.38 AU, is 
located below the RHL of S2, given as 2.05 AU. The same criterion (see Eq. 36) also dictates 
the relatively small width of the radiative GHZ, which extends between 2.16 and 2.66 AU. 
At a binary separation of 2.0 AU, the situation is even more drastic as both the radiative 
CHZ and GHZ are disallowed; the only type of circumbinary habitability remaining is that 
provided by the relatively large EHZ, which is also consistent with the planetary orbital 
stability constraint. 

Additional analyses have been pursued with the binary separation distance 2a taken as 
an independent variable. Hence, for systems of M\ = M2 = 1.0 Mq, it is found that the 
radiative CHZ and GHZ vanish at 2a = 1.52 AU and 2.35 AU, respectively; the associated 
limits of orbital stability are given as 2.78 and 4.31 AU, respectively. Furthermore, the 
distances of fading away for the radiative CHZ and GHZ are identified as 1.76 and 2.45 AU, 
respectively. Therefore, it can be concluded that for equal-mass binary systems with Mi = 
M2 = 1.0 Mq, habitability for relatively wide binaries is lost because of the lack of planetary 
orbital stability occurring at binary separations with the circumbinary RHZ still present. 

We repeated the same type of study for binary systems with Mi = 1.5 Mq and 
M2 = 0.5 Mq. In this case it was found that the radiative CHZ and GHZ vanish at 
2a = 0.54 AU and 1.24 AU, respectively; here the associated limits of orbital stability 
are given at 0.71 and 1.63 AU, respectively. Also, the distances of fading away regarding 
the radiative CHZ and GHZ are found at 1.28 and 2.35 AU, respectively. For this type of 
system it is found that habitability for relatively wide binaries is lost because of the lack 
of existence of the circumbinary RHZ, even though the orbital stability criterion would still 
be met. The fact that the circumbinary RHZ starts to fade away at relatively small binary 
distances, if compared to the case Mi = M2 = 1.0 Mq, is relatively surprising in the view 



of the overall much higher system luminosity because of Mi and L^, oc Mf (e.g., Reid 1987); 
see Sect. 2. However, it is a consequence of the imbalance of binary systems with highly 
unequal distributions of stellar luminosity outputs. 

Radiative imbalance of binary systems in the context of circumbinary habitability is 
strongly related to a property previous discussed, which is that it may occur that the RHL 
for s^^out to be partially or completely located inside of the RHL for s^^n, see Sect. 3.4. 
Therefore, we also report more detailed results about the variation of the inner and outer 
limits (RHLs), which if appropriate amount to the establishment of RHZs; see Eq. (32). 
Specifically, as conveyed through Eq. (36) for a specific RHZ to exist, it is required for the 
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RHL of s^^out to be located completely outside of the RHL for in with the s^^in and the s^^out 
being appropriately paired. 

The condition of Eq. (36) is most readily met in case of equal-mass binary systems with 
small separation distances; e.g., complete separation for the inner and outer limits for the 
appropriate (s^,in, s^.out pairings is found for, e.g.. Mi — M2 — 1.0 Mq and a — 0.25 AU. 
Conversely, the nullification of RHZs is most likely to occur in wide binary systems with 
unequal masses. Figure 7 depicts the ranges of the various RHLs for various cases previously 
discussed. It shows that for Mi = M2 = 1.0 Mq and a = 0.25 AU. the RHLs pertaining to 
the radiative GHZ and GHZ are wiedely separated; thus, the radiative GHZ and GHZ are 
clearly able to exist. 

For Ml = 1.5 Mq, M2 = 0.5 Mq and a = 0.25 AU the pairings {32,84) and (51,55) 
arc again separated, though for (52,54) this is barely the case as the upper limit of the 52 
range almost coincides with the lower limit of the S4 range; nevertheless, the full extents of 
the radiative GHZ and GHZ are permitted, and possibihty of habitabihty for these types of 
RHZs solely depends on the fulfillment of the planetary orbital stabihty constraint. However, 
in case of Mi = 1.5 Mq, M2 = 0.5 Mq and a = 1.0 AU both the radiative GHZ and GHZ 
are disallowed. In this case, only the radiative EHZ is permitted. Fortunately, it is also 
permitted based on the planetary orbital stability constraint (see Fig. 6), implying that this 
type of system offers a tenuous chance of circumbinary habitabihty. 

Although the unequal distribution of stellar mass is identified as highly advantageous 
to the facilitation of planetary orbital stability, in consideration of the lower stability limit 
for P-type orbits for smaller mass ratios /i (see Sect. 4), the situation for the existence of 
the RHZs is much less ideal, even for systems where the stellar primary is highly luminous. 
In this regard the radiative GHZ is in most jeopardy as (54 — 52) constitutes the smallest 
bracket among the various kinds of HZs (see Table 2). More fortunate situations are found 
for the radiative GHZ and EHZ, which are also characterized by larger widths, especially in 
systems of stars with high luminosities. 

5.3. Habitability Classification: Habitabihty Types S, P, ST, and PT 

Another key aspect of this study is provide an appropriate classification of circumstcllar 
habitability, which will be illustrated taking equal-mass binary systems for theoretical main- 
sequence stars of spectral types FO to MO as examples. The required stellar parameters, 
which are the effective temperatures Tefr, the luminosities L*, and the masses of the 
stellar components are given in Table 1. 
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Previously, Dvorak (1982) introduced the terminology of S-type and P-type orbits for 
system planets, which is now widely used by the orbital stability, planetary, and the astrobi- 
ology science communities. Evidently, besides the assessment of orbital stability behaviors, 
these terms are also appropriate for classifying binary system RHZs, if existing. However, 



following previous investigations (e.g.. 


Takeda et al. 


2008 


Dvorak et al.||2010 


et al.||2010 


Kopparapu & Barnes||2010 


Eggl et al.||2012 


), as well as the results 



work, the spatial domain of S-type and P-type habitability depicted by the RHZSs is often 
adversely affected, and in some cases even nullified, by the requirement that system planets 
must be orbitally stable. Thus, if the available extent of the S-type and P-type RHZSs for 
the manifestation of habitability is truncated owing to the additional constraint of planetary 
orbital stability, these zones will be referred to as ST-type and PT-type, respectively, in the 
following. 

Computationally, the occurrence of P, PT, ST, and S-type habitability can be identified 
as follows (see Fig. 8). For sufficiently small binary separations (2a), it is found that both the 
inner and outer limit of the P-type RHZ are located outside of the P-type orbital stability 
limit (see Eq. 43), which constitutes a lower limit of planetary orbital stability (see Sect. 4). 
If the distance of binary separation is increased, the inner and outer limits of the P-type 
RHZ decrease, whereas the P-type orbital stability limit increases. Thus, starting at a certain 
value of 2a, the inner limit of the P-type RHZ will be positioned inside of the P-type orbital 
stability limit; in this case, PT-type habitability is attained. If the binary separation is 
further increased, both the inner and outer limit of the P-type RHZ (if existing) will move 
inside the P-type orbital stability limit; thus, there will be no habitability. 

For mid-sized values of the binary separation distance, the inner and outer limit of 
the S-type RHZ are located inside of the S-type orbital stability limit (see Eq. 42), which 
constitutes an upper limit of planetary orbital stability (see Sect. 4). This setting again entails 
the absence of habitability for binary system planets. If the binary separation distance is 
further increased, the S-type orbital stability will increase and will move beyond the inner 
limit of the S-type RHZ; in this case, ST-type habitability is encountered. Eventually, 
for sufficiently large binary separations, the S-type orbital stability limit is positioned at a 
sufficiently large distance from the stellar component in order for both the inner and outer 
limit of the S-type RHZ, which moved closer to the stellar component, to be located inside 
the S-type orbital stability limit. Now, S-type habitability is found because the total width 
of the S-type RHZ is available for facilitating habitability. 

S-type, P-type, ST-type, and PT-type habitability for the CHZ, GHZ, and EHZ for 
equal-mass main-sequence star binaries of spectral types FO to MO are displayed in Tables 8, 
9, and 10. The CHZ, GHZ, and EHZ correspond to the habitability limits (^2,54), (si,S5), 
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and {si,Se), respectively, see Sect. 3.4. It is found that for the evaluated array of spectral 
type / binary separation combinations, with the latter inspected between distances (2a) of 
0.2 and 10.0 AU, the CHZ exhibits the most cases of non-existent habitability, as expected, 
owing to its property that it has the smallest extent compared to the GHZ and EHZ. But if 
habitability for the CHZ is encountered, it is usually of type S or P rather than ST or PT, 
noting that because of the CHZ's small extent, cutting through (i.e., truncation) of the RHZ 
by the S-type or P-type orbital stabihty limit (see Eqs. 42 and 43) is relatively rare. 

On the other hand, for the GHZ and the EHZ, owing to their significantly larger ex- 
tents, there is a much smaller range of binary separation distances without the possibility of 
circumstellar habitability. However, if habitability is realized, there is a higher percentage of 
PT versus P cases, as well as ST versus S cases, found for the GHZ, and even more so for the 
EHZ, compared to the CHZ. This behavior is again due to the fact that the GHZ, and even 
more so the EHZ, is much more extended than the EHZ. Thus, there is a higher likelihood 
that the S-typc or P-type planetary orbital stability limit interferes with those zones. Most 
importantly, there is also no case of an equal-star binary system between the spectral types 
FO and MO where both S/ST or P/PT-type habitability is found to be possible. This result 
has already been suggested by the Eqs. (14) and (16), see Sect. 3.2. 

For G2 V stars, i.e., two Suns in a circular orbit, the following picture emerges: For the 
CHZ, P-type habitability is encountered up to a distance of 0.6 AU, followed by PT-type 
habitability up to a distance (2a) of 0.8 AU. ST-type habitability is encountered at a binary 
separation of 8.3 AU, whereas S-typc habitability occurs at 12.0 AU. Regarding the GHZ 
and EHZ, somewhat smaller values compared to the CHZ occur for P-type habitability, and 
moderately larger values occur for PT-type habitability. The biggest difference between 
the GHZ and EHZ, on one hand, and the CHZ, on the other hand, is found for S-type 
habitabihty: For the GHZ and EHZ, this mode of habitability requires binary separations 
of 14.6 and 21.0 AU, respectively; see Table 11 for complete results. 

Another representation of our results is given in Fig. 9. It conveys a depiction of S/ST- 
type and P/PT-type habitability using different gray scales regarding the CHZ, GHZ, and 
EHZ. It also shows a notable domain with no habitability at all. Note that no separate 
S/ST- type contour is attained for the EHZ, owing to the fact that the inner boundary of 
the EHZ (i.e., si; see Table 2) agrees with the inner boundary of the GHZ, thus rendering 
the same criterion for the inner limit of habitability. For the distance limits (2a) of ST-type 
habitability at the high-luminosity end of the stellar range (attained outside the figure's 
frame), the values for the F5 V, F2 V, and FO V stars pertaining to the GHZ and EHZ are 
given as 11.7, 13.6, and 15.9 AU, whereas for the CHZ, they are given as 13.3, 15.5, and 
18.1 AU, respectively. 
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6. Summary and Conclusions 



This study is aimed at developing a new method for a comprehensive assessment of S- 
type and P-type habitabihty in stellar binary systems. P-type orbits occur when the planet 
orbits both binary components, whereas in case of S-type orbits the planet orbits only one of 
the binary components with the second component considered a perturbator. An important 
characteristic of the method is that the joint constraint of orbital stability for a system planet 
as well as the necessity for providing a habitable environment for that planet through the 
stellar radiative energy fluxes ("radiative habitable zone") needs to be met, a property also 



recognized in previous studies (e.g., Takeda et al. 2008 Kopparapu & Barnes 2010 Eggl et 
iL] [20l2| [Kane fc Hinkd|[20l3| ) . 



In conjunction with the present study, a habitability classification of type S, P, ST, and 
PT is also proposed, which considers whether or not S-type and P-type radiative habitable 
zones are truncated based on the additional constraint of planetary orbital stability. This 
classification scheme may be of general relevance as it can be applied to both equal-mass 
and non-equal mass binary systems, as well as to systems with binaries in elliptical orbits. 

A key aspect of the proposed method is the introduction of a combined algebraic formal- 
ism for the assessment of both S-type and P-type habitability; in particular, mathematical 
criteria are presented allowing to determine for which systems S-type and P-type RHZs are 
realized. In this regard, apriori choices about the presence of S-type and P-type RHZs are 
neither necessary nor possible as the existence of S-type as well as P-type RHZs is prolif- 
erated through well-defined mathematical conditions regarding the underlying descriptive 
fourth-order algebraic equation. The coefficients of the polynomial are given by the binary 
separation distance (2a), the solar system-based parameter for the limit of habitability (se), 
and the modified values for the luminosities {L[, L'2) of the stellar binary components. Re- 
garding the binary system habitable zone, we consider conservative, general and extended 
zones of habitability, noting that their inner and outer limits are informed by previous solar 



system investigations (e.g., Kasting et al. 1993; Underwood et al. 2003; Selsis et al. 2007). 



After combining the requirements pertaining to the existence of the RHZs with the ne- 
cessity of orbital stability for system planets, we were able to identify five different cases, 
which are: S-type and P-type habitability provided by the full extent of the RHZ; habitabil- 
ity, where the RHZ is truncated by the additional constraint of planetary orbital stability 
(labelled as ST and PT-type, respectively); and cases of no habitability at all. This will allow 
a general classification of habitability which is expected to be also useful for higher-order 
systems, i.e., N > 2. We provided examples for various cases, i.e., equal- mass and nonequal- 
mass systems, by focusing on theoretical main-sequence stars based on parameters given by 



Castelli &; Kurucz (2004), Gray (2005), and Kurucz (2005). As a limitation of scope, we 
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restricted our investigation to binary systems with the stellar components in circular orbits; 
note that systems in elliptical orbits will be examined in Paper II. 

In the segment of applications, we examined the existence of habitable S-type orbits 
for selected examples. We found that regarding the RHZs, owing to the relatively large 
separation of the stellar components, the effect of the stellar secondary on the extents of the 
RHZs is typically very minor. The secondary's main influence on circumstellar habitability 
consists in imposing restrictions regarding planetary orbital stability in form of an upper 
stability limit around each stellar component, which often truncates or nullifies S-type plan- 
etary habitability. Previous studies, which also include binary systems in elliptical orbits. 



thus entailing much more intricate outcomes, were given by, e.g., Eggl et al. (2012). In the 
framework of our study, we specifically considered the radiative EHZ. It was found that this 
kind of zone is most affected by the limitation of planetary orbital stability as it is closest 
to the secondary stellar component. 

More complicated scenarios emerge for potentially habitable P-type orbits, which were 
studied in detail through comparing pairs with masses of (1.0 Mq, 1.0 Mq) to pairs with 
masses of (1.5 Mq, 0.5 Mq). Results for pairs of (1.25 Mq, 0.75 Mq) were given as in- 
termediate case for additional information. The overall focus of this effort was to deduce 
information on how the distribution of stellar mass and, by implication, stellar luminosity 
impacts the habitability prospects of main-sequence star binary systems. It was found that 
although the unequal distribution of stellar mass is highly advantageous to the facilitation 
of planetary orbital stability, in consideration of the lower stability limit for P-type orbits 
(see Sect. 4), the situation for the existence of the RHZs may be highly complicated. 

The reason is that an unequal distribution of mass always entails un unequal distribu- 
tion of luminosity with the more massive stellar component being overwhelmingly luminous 
in accord to the mass-luminosity relationship, i.e., oc Mf ( |Rddl[l987| ; see Sect. 2. Gener- 



ally, this property drastically increases the overall habitability prospects of a P-type binary 
system. However, there is an important complication that is particularly relevant for rel- 
atively wide binaries: Careful evaluations show that for some systems the outer limit of 
distinct RHZs may be positioned inside the inner limit of that zone, at least in some parts of 
the system (see, e.g.. Fig. 8) , which effectively nullifies that zone, irrespectively of planetary 
orbital stability considerations. In this regard, the GHZ is in highest jeopardy as (54 — S2) 
constitutes the smallest bracket among the various kinds of HZs (see Table 2). 

Future work will deal with a significant augmentation of our method to other systems, 
including systems with binary components in elliptical orbits (see Paper II). This will allow 
us to compare the potential and outcome of our method, as well as results for individual 
systems, to other findings in the literature. We also expect our method to be applicable to 



- 28 - 



general binary systems with main-sequence stars as well as to systems containing evolved 
stars; this latter effort is motivated by observational evidence and supporting theoretical 
efforts indicating that planets are also able to exist around stars that have left the main- 
sequence (e.g., |Sato et~aL]|2003| |Ramm et all2009[ [Eberle fc Cuntz||2010[ [Doyle et al.||201lD . 
Particularly, it is highly desirable to augment our method to systems of higher order, as 
motivated by the steady progress in theory as well as ongoing and future observational 
discoveries of exosolar planetary systems. 
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habitability in binary systems based on the developed method. 



- 29 - 



REFERENCES 

Abramowitz, M., & Stegun, I. A. 1972, Handbook of Mathematical Functions with Formulas, 
Graphs, and Mathematical Tables, 9th printing (New York: Dover) 

Asghari, N., Broeg, C, Carone, L., et al., 2004, A&A, 426, 353 

Beyer, W. H. 1987, Handbook of Mathematical Sciences, 6th edn. (Boca Raton: CRC Press) 
Bonavita, M., & Desidera, S. 2007, A&A, 468, 721 

Bronshtein, I. N., & Semendyayev, K. A. 1997, Handbook of Mathematics, 3rd edn. (New 
York: Springer) 

Castelh, F., & Kurucz, R. L. 2004, in lAU Symp. 210, Modelling of Stellar Atmospheres ed. 
N. E. Piskunov, W. W. Weiss, & D. F. Gray (San Francisco: ASP), CD-ROM Poster 
20; [a'rXiv:astro-ph/0405087| 

Cuntz, M., & Yeager, K. E. 2009, ApJ, 697, L86 

Cuntz, M., Eberle, J., & Musielak, Z. E. 2007, ApJ, 669, L105 

Doyle, L. R., Carter, J. A., Fabrycky, D. C, et al. 2011, Science, 333, 1602 

Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 

Dvorak, R. 1982, Osterreichische Akademie der Wissenschaften, Mathematisch- 
Naturwissenschafthche Klasse, Sitzungsberichte Abt. 2, 191 (10), 423 

Dvorak, R. 1984, Celest. Mech. Dyn. Astron., 34, 369 

Dvorak, R. 1986, A&A, 167, 379 

Dvorak, R., Pilat-Lohinger, E., Bois, E., et al. 2010, Astrobiology, 10, 33 
Eberle, J., & Cuntz, M. 2010, ApJ, 721, L168 
Eberle, J., Cuntz, M., & Musielak, Z. E. 2008, A&A, 489, 1329 
Eggenberger, A., Udry, S., & Mayor, M. 2004, A&A, 417, 353 

Eggenberger, A., Udry, S., Chauvin, G., Beuzit, J.-L., Lagrange, A.-M., Segransan, D., & 
Mayor, M. 2007, A&A, 474, 273 

Eggl, S., Pilat-Lohinger, E., Georgakarakos, N., Gyergyovits, M., & Funk, B. 2012, ApJ, 
752, A74 



- 30 - 

Fatuzzo, M., Adams, F. C, Gauvin, R., & Proszkow, E. M. 2006, PASP, 118, 1510 
Forget, F., & Pierrehumbert, R. T. 1997, Science, 278, 1273 

Gray, D. F. 2005, The Observation and Analysis of Stellar Photospheres, 2nd edn. (Cam- 
bridge: Cambridge University Press) 

Haghighipour, N., Dvorak, R., & Pilat-Lohinger, E. 2010, in Planets in Binary Systems, 
ed. N. Haghighipour, Astrophysics and Space Science Library, Vol. 366 (New York: 
Springer Science + Business Media), p. 285 

Halevy, I., Pierrehumbert, R. T., & Schrag, D. P. 2009, J. Geophys. Res., 114, D18112 

Holman, M. J., & Wiegert, P. A. 1999, AJ, 117, 621 

Jackson, J. D. 1999, Classical Electrodynamics, 3rd edn., John Wiley & Sons, Inc., Hoboken, 
NJ 

Jones, B. W., Sleep, P. N., & Chambers, J. E. 2001, A&A, 366, 254 
Jones, B. W., & Sleep, P. N. 2010, MNRAS, 407, 1259 
Kane, S. R., & Hinkel, N. R. 2013, ApJ, 762, A7 

Kasting, J. F., Whitmire, D. P., & Reynolds, R. T. 1993, Icarus, 101, 108 

Kopparapu, R. K., & Barnes, R. 2010, ApJ, 716, 1336 

Kurucz, R. L. 2005, Mem. S. A. It., 8, 14 

Lada, C. J. 2006, ApJ, 640, L63 

Menou, K., & Tabachnik, S. 2003, ApJ, 583, 473 

Mischna, M. A., Kasting, J. F., Pavlov, A., & Freedman, R. 2000, Icarus, 145, 546 
Musielak, Z. E., Cuntz, M., Marshall, E. A., & Stuit, T. D. 2005, A&A, 434, 355 
Noble, M., Musielak, Z. E., & Cuntz, M. 2002, ApJ, 572, 1024 
Orosz, J. A., Welsh, W. F., Carter, J. A., et al. 2012, Science, 337, 1511 
Patience, J., White, R. J., Chez, A. M., et al. 2002, ApJ, 581, 654 

Ferryman, M. 2011, The Exoplanet Handbook (Cambridge: Cambridge University Press) 



Quarles, B., Eberle, J., Musielak, Z. E., & Cuntz, M. 2011, A&A, 533, A2 
Quarles, B., Musielak, Z. E., & Cuntz, M. 2012, ApJ, 750, A14 

Quintana, E. V., & Lissauer, J. J. 2010, in Planets in Binary Systems, ed. N. Haghighipour, 
Astrophysics and Space Science Library, Vol. 366 (New York: Springer Science + 
Business Media), p. 265 

Raghavan, D., Henry, T. J., Mason, B. D., et al. 2006, ApJ, 646, 523 

Raghavan, D., McAlister, H. A., Henry, T. J., et al. 2010, ApJS, 190, 1 

Ramm, D. J., Pourbaix, D., Hearnshaw, J. B., & Komonjinda, S., 2009, MNRAS, 394, 1695 

Reid, N. 1987, MNRAS, 225, 873 

Sandor, Zs., Suli, A., Erdi, B., Pilat-Lohinger, E., Dvorak, R. 2007, MNRAS, 375, 1495 
Sato, B., et al., 2003, ApJ, 597, L157 

Satyal, S., Quarles, B., & Hinse, T. C. 2013, MNRAS, submitted 

Selsis, P., Kasting, J. P., Levrard, B., Paillet, J., Ribas, I., & Delfosse, X. 2007, A&A, 476, 
1373 

Setiawan, J., Hatzes, A. P., von der Liilie, et al. 2003, A&A, 398, L19 

Slawson, R. W., Prsa, A., Welsh, W. P., et al. 2011, AJ, 142, 160 

Stix, M. 2004, The Sun: An Introduction, 2nd edn. (A&A Library, Berhn: Springer) 

Szenkovits, P., & Mako, Z. 2008, Celest. Mech. Dyn. Astron., 101, 273 

Takeda, G., Kita, R., & Rasio, P. A. 2008, ApJ, 683, 1063 

Trilling, D. E., Stansberry, J. A., Stapelfeldt, K. R., et al. 2007, ApJ, 658, 1289 

Underwood, D. R., Jones, B. W., & Sleep, P. N. 2003, Int. J. AstrobioL, 2, 289 

Wilhams, D. M., & Pollard, D. 2002, Int. J. AstrobioL, 1, 61 



This preprint was prepared with the AAS IM^jX macros v5.0. 



- 32 - 





Fig. 1. — Set-up for the mathematical treatment of S-type (top) and P-type (bottom) hab- 
itable zones of binary systems as given by the stellar radiative fluxes. Note that the stars SI 
and S2 have been depicted as identical for convenience. 
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S-Type and P-Type Habitability for Equal Binaries 




Distance (AU) 

Fig. 2. — Habitability limits for equal- mass binaries regarding S-type and P-type planetary 
orbits pertaining to the RHZ. The various limits, separated by grayish areas, refer to different 
values of given as 0.84, 0.95, 1.00, 1.37, 1.67, and 2.40 AU (from left to right; corresponding 
to £ = 1 to 6); see Sect. 2 and 3.2 for details. 





Fig. 3. — Examples of S-type and P-type RHZs for different systems obtained through 
solutions of the governing fourth-order polynomial equations. The thick solid lines indicate 
the RHLs, corresponding to the inner and outer limit of habitability based on = 0.84 
and 1.67 AU, respectively. Note that in the S-type case the RHLs are bended toward the 
center, whereas in the P-type case they are of elliptical shape. The top row displays systems 
with L'^ = L'2 = Lq, whereas the bottom top row displays systems with L[ = 1.5 Lq and 
L2 = 0.5 Lq. The left column assumes separation distances (2a) of 5.0 AU rendering P-type 
RHZs, whereas the right column assumes separation distances of 0.5 AU resultant in S-type 
RHZs. The gray areas indicate the appropriate circular regions (annuli), referred to as RHZs, 
for each case. 
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Maximum Distance for P-Type Habitability 
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Recast Luminosity L'(L0) 



Fig. 4. — Maximum of permissible separation distance 2a of equal-mass binary components, 
which still permits the possibility of P-type RHZs (see Eq. 41). Results are given for the 
CHZ, GHZ, and EHZ. 
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Fig. 5. — Domains of S-type habitability for main-sequence star binary systems. The left 
column depicts systems with masses of Mi = M2 = 1.0 Mq, whereas the right column 
depicts systems with masses of Mi = 1.5 Mq and M2 = 1.0 Mq. The separation distances 
(2a) are given as 10 AU (top) and 20 AU (bottom). Results for the conservative, general 
and extended RHZs are given as dark gray, medium gray, and light gray areas, respectively. 
The orbital stability hmits are indicated by red dashed lines. 
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Fig. 6. — Domains of P-type habitability for main-sequence star binary systems. The left 
column depicts systems with masses of Mi = M2 = 1.0 Mq, whereas the right column 
depicts systems with masses of Mi = 1.5 Mq and M2 = 0.5 Mq. The separation distances 
(2a) are given as 0.5 AU (top), 1.0 AU (middle), and 2.0 AU (bottom). Results for the 
conservative, general and extended RHZs are given as dark gray, medium gray, and light 
gray areas, respectively. The orbital stability limits are indicated by red dashed lines. 
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a = 0.25 AU M =1.0M M=1.0M 

1 o 2 o 



a = 0.25 AU M =1.5M M=0.5M 

1 o 2 o 





a = 1.0AU M =1.0M M=1.0M 

1 o 2 o 



a = 1 .0 AU M =1 .5M M =0.5M 

1 o 2 o 





and 



Fig. 7. — Ranges of the inner and outer limits of RHZs, given as Mm(7l{z,ip) 

Max(n{z,ip)) , pertaining to se = 0.84, 0.95, 1.37, 1.67, and 2.40 AU (from left to right; 

see labels 1, 2, 4, 5, and 6). The figure panel depicts binary systems with separation distances 
(2a) of 0.5 AU (top) and 2.0 AU (bottom) and with masses of Mi = M2 = 1.0 Mq (left) and 
Ml = 1.5 Mq and M2 = 0.5 Mq (right); see also Tables 5 to 7 for additional information. 
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6 9 
Distance (AU) 

Fig. 8. — Facilitation of habitability in an equal-star binary system of two G2 V stars as 
function of the binary separation distance (2a) pertaining to the GHZ, i.e., s^^in = 0.84 and 
se,out = 1-67 AU. The two red lines indicate the limits of the P-type RHZ (i.e., RHZin and 
RHZout; see Eqs. (40a) and (40b)), whereas the two blue lines indicate the limits of the 
S-type RHZ (see Eqs. (37a) and (37b)). The available S-type and P-type RHZs are depicted 
as grayish areas. Furthermore, the cyan dashed line indicates the P-type orbital stability 
limit, whereas the violet dashed line indicates the S-type orbital stability limit; see Eqs. (43) 
and (42), respectively. Note that the P-type orbital stability limit constitutes a lower limit 
(i.e., planetary orbital stability is attained in the domain to the upper left of that line), 
whereas the S-type orbital stability limit constitutes an upper limit (i.e., planetary orbital 
stability is attained below that line). Therefore, the intersections between the P-type RHZs 
and the P-type orbital stability limit allow to identify the domains of P-type and PT-type 
habitability given to the left of the intersection points. Also, the intersections between the 
S-type RHZs and the S-type orbital stability limit allow to identify the domains of S-type 
and ST-type habitability given to the right of the intersection points. No habitability is 
found between the intersection points PT and ST. Similar figures hold for other types of 
binary systems, including systems of non-equal stellar masses. 
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FO F5 GO G5 KO K5 MO 

Spectral Type 



Fig. 9. — Depiction of S/ST-type and P/PT-type habitability based on the joint constraint 
of planetary orbital stability and the availability of a habitable environment provided by the 
stellar radiative energy fluxes. The domains of the conservative, general and extended HZs 
are depicted as dark gray, medium gray, and light gray areas, respectively. The results refer 
to the particular case of equal-mass binaries in circular orbits, depicted as main-sequence 
stars of identical spectral types. 
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Table 1. Stellar Parameters 



Sp. Type 








'S'rel,* 






(K) 


(Re) 


(Le) 




(Me) 


FO 


7178 


1.62 


6.255 


1.145 


1.60 


F2 


6909 


1.48 


4.481 


1.113 


1.52 


F5 


6528 


1.40 


3.196 


1.072 


1.40 


F8 


6160 


1.20 


1.862 


1.037 


1.19 


GO 


5943 


1.12 


1.405 


1.019 


1.05 


G2 


5811 


1.08 


1.194 


1.009 


0.99 


G5 


5657 


0.95 


0.830 


0.997 


0.91 


G8 


5486 


0.91 


0.673 


0.985 


0.84 


KO 


5282 


0.83 


0.481 


0.971 


0.79 


K2 


5055 


0.75 


0.330 


0.957 


0.74 


K5 


4487 


0.64 


0.149 


0.926 


0.67 


K8 


4006 


0.53 


0.066 


0.905 


0.58 


MO 


3850 


0.48 


0.045 


0.900 


0.51 



Note, calculated using se ^ S3. 
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Table 2. Habitability Limits for Solar Environments 





Si 

(AU) 


Description 


1 


0.84 


Runaway greenhouse effect 


2 


0.95 


Start of water loss 


3 


1.00 


Earth-equivalent position 


4 


1.37 


First CO2 condensation 


5 


1.67 


Maximum greenhouse effect, no clouds 


6 


2.40 


Maximum greenhouse effect, 100% clouds 



Note. — See text for references. 



Table 3. Habitable Zones of Single Main-Sequence Stars 



Sp. Type 




Habitable Zone 






HZ-iG 


HZ-iC 


HZ-oC 


HZ-oG 


HZ-oE 




(AU) 


(AU) 


(AU) 


(AU) 


(AU) 


FO 


1.978 


2.253 


2.909 


3.659 


5.485 


F2 


1.696 


1.929 


2.543 


3.178 


4.723 


F5 


1.456 


1.653 


2.243 


2.780 


4.085 


F8 


1.128 


1.278 


1.782 


2.192 


3.188 


GO 


0.987 


1.118 


1.584 


1.940 


2.805 


G2 


0.915 


1.035 


1.480 


1.808 


2.606 


G5 


0.766 


0.866 


1.253 


1.527 


2.192 


G8 


0.694 


0.784 


1.148 


1.395 


1.994 


KO 


0.590 


0.667 


0.990 


1.199 


1.705 


K2 


0.492 


0.555 


0.837 


1.010 


1.429 


K5 


0.335 


0.378 


0.593 


0.709 


0.990 


K8 


0.224 


0.253 


0.410 


0.487 


0.674 


MO 


0.187 


0.211 


0.345 


0.409 


0.565 



Note. — HZ-i and HZ-o indicate the inner and outer 
limits, respectively, of single main-sequence star habitable 
zones. The added letters C, G and E are in reference to 
the conservative, general and extended habitable zones, 
respectively (see text for details). 
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Table 4. Target List, Sections 5.1 and 5.2 





TeS 




L' 


HZ^ eqv 


(Mo) 


(K) 


(Le) 


(Le) 


(AU) 


0.50 


3857 


0.045 


0.050 


0.22 


0.75 


5100 


0.357 


0.372 


0.61 


1.00 


5833 


1.228 


1.216 


1.10 


1.25 


6256 


2.173 


2.078 


1.44 


1.50 


6842 


4.233 


3.830 


1.96 



Note. — The position of HZ®^, 
corresponds to — s^. 
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Table 5. P-Type Stellar Habitability; a = 0.25 AU 



HZ 



Binary System Data 



(AU) 



(AU) 



(AU) 



(AU) 



Model: 


Ml = 1.0 Mq, M2 


= 1.0 Mq 




0.84 


1.29 


1.38 


0.92 


0.95 


1.46 


1.54 


0.92 


1.00 


1.54 


1.62 


0.92 


1.37 


2.10 


2.16 


0.92 


1.67 


2.58 


2.62 


0.92 


2.40 


3.72 


3.76 


0.92 



Model: iVh = 1-25 Mq, M2 = 0.75 M( 



0.84 1.19 1.52 0.79 

0.95 1.36 1.69 0.79 

1.00 1.45 1.79 0.79 

1.37 1.96 2.28 0.79 

1.67 2.43 2.76 0.79 

2.40 3.58 3.91 0.79 



Model: 


Ml = 1.5 Mq, M2 


= 0.5 Mq 




0.84 


1.42 


1.91 


0.66 


0.95 


1.65 


2.14 


0.66 


1.00 


1.73 


2.22 


0.66 


1.37 


2.28 


2.76 


0.66 


1.67 


2.90 


3.38 


0.66 


2.40 


4.40 


4.89 


0.66 
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Table 6. P-Type Stellar Habitability; a = 0.50 AU 



HZ 



Binary System Data 



(AU) 



(AU) 



(AU) 



(AU) 



Model: 


Ml = 1.0 Mq, M2 


= 1.0 Mq 




0.84 


1.21 


1.54 


1.83 


0.95 


1.40 


1.69 


1.83 


1.00 


1.48 


1.76 


1.83 


1.37 


2.06 


2.28 


1.83 


1.67 


2.54 


2.72 


1.83 


2.40 


3.70 


3.83 


1.83 



Model: iVh = 1-25 Mq, M2 = 0.75 M( 



0.84 1.11 1.75 1.57 

0.95 1.29 1.92 1.57 

1.00 1.38 2.02 1.57 

1.37 1.90 2.49 1.57 

1.67 2.35 2.96 1.57 

2.40 3.46 4.12 1.57 



Model: 


Ml = 1.5 Mq, M2 


= 0.5 Mq 




0.84 


1.21 


2.16 


1.32 


0.95 


1.43 


2.38 


1.32 


1.00 


1.51 


2.47 


1.32 


1.37 


2.05 


3.00 


1.32 


1.67 


2.66 


3.62 


1.32 


2.40 


4.17 


5.13 


1.32 
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Table 7. P-Type Stellar Habitability; a = 1.0 AU 



HZ 



Binary System Data 



(AU) 



(AU) 



(AU) 



(AU) 



Model: 


Ml = 1.0 Mq, M2 


= 1.0 Mq 




0.84 


0.85 


1.98 


3.66 


0.95 


1.10 


2.11 


3.66 


1.00 


1.20 


2.18 


3.66 


1.37 


1.87 


2.64 


3.66 


1.67 


2.39 


3.05 


3.66 


2.40 


3.60 


4.09 


3.66 



Model: Mi = 1.25 Mq, M2 = 0.75 M( 



0.84 0.70 2.23 3.15 

0.95 0.95 2.40 3.15 

1.00 1.07 2.50 3.15 

1.37 1.69 2.95 3.15 

1.67 2.18 3.42 3.15 

2.40 3.32 4.55 3.15 



Model: 


Ml = 1.5 Mq, M2 


= 0.5 Mq 




0.84 


0.83 


2.66 


2.63 


0.95 


1.09 


2.88 


2.63 


1.00 


1.18 


2.96 


2.63 


1.37 


1.74 


3.50 


2.63 


1.67 


2.28 


4.12 


2.63 


2.40 


3.71 


5.63 


2.63 
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Table 8. Assessment of the Conservative Habitable Zone of Equal-Mass Binary Systems 



bp. iype 






Binary Separation Distance (2a) 








0.2 


0.5 


1.0 


2.0 3.0 4.0 


5.0 


7.5 


10.0 




(AU) 


(AU) 


(AU) 


(AU) (AU) (AU) 


(AU) 


(AU) 


(AU) 


FO 


P 


P 


p 

r 










Fl 


P 


P 


r 










F2 


P 


P 


p 
r 










F3 


P 


P 


p 
r 










F4 


P 


P 


p 










F5 


P 


P 


p 

r 










F6 


P 


P 


r 1 










F7 


P 


P 


r 1 










F8 


P 


P 


PT 










F9 


P 


P 










ST 


GO 


P 


P 










ST 


Gl 


P 


P 










ST 


G2 


P 


P 










ST 


G3 


P 


P 










ST 


G4 


P 


P 








ST 


ST 


G5 


P 


P 








ST 


ST 


G6 


P 


P 








ST 


S 


G7 


P 


P 








ST 


s 


G8 


P 


PT 








ST 


s 


G9 


P 


PT 








ST 


s 


KO 


P 


PT 








ST 


s 


Kl 


P 


PT 






ST 


S 


s 


K2 


P 








ST 


s 


s 


K3 


P 








ST 


s 


s 


K4 


P 






ST 


ST 


s 


s 


K5 


P 






ST 


S 


s 


s 


K6 


P 






ST ST 


s 


s 


s 


K7 


PT 






ST S 


s 


s 


s 


K8 


PT 






ST S 


s 


s 


s 
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Table 8 — Continued 



Sp. Type 






Binary Separation Distance (2a) 








0.2 


0.5 


1.0 2.0 3.0 4.0 5.0 


7.5 


10.0 




(AU) 


(AU) 


(AU) (AU) (AU) (AU) (AU) 


(AU) 


(AU) 


K9 


PT 




ST ST S S 


S 


s 


MO 


PT 




ST S S S 


s 


s 



Note. — S and P indicate S-type and P-typc habitability of the binary system, 
respectively, given by the full spatial extent of the radiative habitable zones. ST 
and PT indicate S-type and P-type habitability of the binary system, respectively, 
that is truncated by the additional constraint of orbital stability for a putative 
system planet. 
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Table 9. Assessment of the General Habitable Zone of Equal-Mass Binary Systems 



bp. iype 






Binary Separation Distance (2a) 








0.2 


0.5 


1.0 


2.0 


3.0 4.0 


5.0 


7.5 


10.0 




(AU) 


(AU) 


(AU) 


(AU) 


(AU) (AU) 


(AU) 


(AU) 


(AU) 


FO 


P 


P 


p 

r 


PT 










Fl 


P 


P 


r 












F2 


P 


P 


p 
r 












F3 


P 


P 


p 
r 












F4 


P 


P 


PT 
r ± 












F5 


P 


P 


r 1 












F6 


P 


P 


PT 












F7 


P 


P 


PT 










ST 


F8 


P 


P 


PT 










ST 


F9 


P 


P 


PT 










ST 


GO 


P 


P 


PT 










ST 


Gl 


P 


P 


PT 










ST 


G2 


P 


P 


PT 








ST 


ST 


G3 


P 


P 










ST 


ST 


G4 


P 


P 










ST 


ST 


G5 


P 


PT 










ST 


ST 


G6 


P 


PT 










ST 


ST 


G7 


P 


PT 










ST 


ST 


G8 


P 


PT 










ST 


ST 


G9 


P 


PT 










ST 


ST 


KO 


P 


PT 








ST 


ST 


S 


Kl 


P 


PT 








ST 


ST 


s 


K2 


P 


PT 






ST 


ST 


ST 


s 


K3 


P 


PT 






ST 


ST 


ST 


s 


K4 


P 








ST 


ST 


S 


s 


K5 


P 








ST ST 


ST 


s 


s 


K6 


PT 








ST ST 


S 


s 


s 


K7 


PT 








ST ST 


s 


s 


s 


K8 


PT 






ST 


ST S 


s 


s 


s 
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Table 9 — Continued 



Sp. Type 






Binary Separation Distance (2a) 








0.2 


0.5 


1.0 2.0 3.0 4.0 5.0 


7.5 


10.0 




(AU) 


(AU) 


(AU) (AU) (AU) (AU) (AU) 


(AU) 


(AU) 


K9 


PT 




ST ST S S 


S 


s 


MO 


PT 




ST ST S S 


s 


s 



Note. — 



See Table 8 for comments. 
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Table 10. Assessment of the Extended Habitable Zone of Equal-Mass Binary Systems 



bp. iype 






Binary Separation Distance (2a) 








0.2 


0.5 


1.0 


2.0 


3.0 


4.0 


5.0 


7.5 


10.0 




(AU) 


(AU) 


(AU) 


(AU) 


(AU) 


(AU) 


(AU) 


(AU) 


(AU) 


FO 


P 


P 


p 

r 


PT 


PT 










Fl 


P 


P 


r 


PT 












F2 


P 


P 


p 
r 


PT 












F3 


P 


P 


p 
r 


PT 












F4 


P 


P 


PT 
r ± 


PT 












F5 


P 


P 


r 1 


PT 












F6 


P 


P 


PT 


PT 












F7 


P 


P 


PT 


PT 










ST 


F8 


P 


P 


PT 












ST 


F9 


P 


P 


PT 












ST 


GO 


P 


P 


PT 












ST 


Gl 


P 


P 


PT 












ST 


G2 


P 


P 


PT 
r i 










ST 


ST 


G3 


P 


P 


PT 










ST 


ST 


G4 


P 


P 


PT 










ST 


ST 


G5 


P 


PT 


PT 










ST 


ST 


G6 


P 


PT 


PT 










ST 


ST 


G7 


P 


PT 


PT 










ST 


ST 


G8 


P 


PT 


PT 










ST 


ST 


G9 


P 


PT 


PT 










ST 


ST 


KO 


P 


PT 










ST 


ST 


ST 


Kl 


P 


PT 










ST 


ST 


ST 


K2 


P 


PT 








ST 


ST 


ST 


ST 


K3 


P 


PT 








ST 


ST 


ST 


ST 


K4 


P 


PT 








ST 


ST 


ST 


S 


K5 


P 


PT 






ST 


ST 


ST 


ST 


s 


K6 


PT 








ST 


ST 


ST 


S 


s 


K7 


PT 








ST 


ST 


ST 


s 


s 


K8 


PT 






ST 


ST 


ST 


ST 


s 


s 
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Table 10— Continued 



Sp. Type 






Binary Separation Distance (2a) 








0.2 


0.5 


1.0 2.0 3.0 4.0 5.0 


7.5 


10.0 




(AU) 


(AU) 


(AU) (AU) (AU) (AU) (AU) 


(AU) 


(AU) 


K9 


PT 




ST ST ST S 


S 


s 


MO 


PT 




ST ST ST S 


s 


s 



Note. — 



See Table 8 for comments. 
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Table 11. Limits of habitability in a G2 V equal-star binary system 



HZ Type 


P 


PT 


ST 


S 




(AU) 


(AU) 


(AU) 


(AU) 


GHZ 


0.6 


0.8 


8.3 


12.0 


GHZ 


0.5 


1.0 


7.4 


14.6 


EHZ 


0.5 


1.5 


7.4 


21.0 



Note. — The limits of P and PT consti- 
tute upper limits, whereas limits of S and 
ST constitute lower limits. No habitability 
is found between the limits PT and ST. 



